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SUMMARY 

A mixed finite-difference scheme is presented for the free -vibration analysis of 
simply supported closed noncircular cylindrical shells. The problem is formulated in 
terms of eight first-order differential equations in the circumferential coordinate which 
possess a symmetric coefficient matrix and are free of the derivatives of the elastic and 
geometric characteristics of the shell. In the finite-difference discretization, two inter- 
lacing grids are used for the different fundamental unknowns in such a way as to avoid 
averaging in the difference-quotient expressions used for the first derivative. The 
resulting finite -difference equations are symmetric. The inverse-power method is used 
for obtaining the eigenvalues and eigenvectors. 

Numerical studies of the effects of reducing the local discretization error and of 
mesh refinement on the accuracy and convergence of the solutions obtained by the scheme 
developed, as well as by the conventional schemes, are discussed. Both oval and elliptic 
profiles with constant and variable thicknesses have been considered, and, in all cases, 
a monotonic convergence for the eigenvalues was obtained with the reduction of the local 
discretization error and/or the increase in the number of finite -difference intervals in 
the case of the modified scheme. Comparisons were also made between the results 
obtained from this study and results from some of the previous approximate analyses. 

It was found that the proposed scheme, in addition to a number of other advantages, leads 
to highly accurate results even when a small number of finite-difference intervals are 
used. 


INTRODUCTION 


Although considerable literature has been devoted to the free -vibration analysis of 
circular cylindrical shells, investigations of the vibrations of noncircular cylindrical 
shells are rather limited in extent. Cylindrical shells having a noncircular cross section 
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are used in deep submersibles, aerospace vehicles, and other industrial applications, 
and, therefore, an understanding of their vibration characteristics is desirable. Because 
the radius of curvature of the noncircular cylinder varies with the circumferential coor- 
dinate, closed-form or analytic solutions cannot be obtained, in general, for this class of 
shells. Numerical or approximate techniques are necessary for their analysis. 

Various approximate methods have been used for the free -vibration analysis of 
noncircular cylinders. Perturbation techniques were used by Klosner (refs. 1 and 2) 
for infinitely long cylinders and by Malkina (ref. 3) for shells with simply supported 
curved edges. Ivanyuta and Finkel'shteyn (ref. 4) applied the variational method of 
Oniashvili (ref. 5) in their study of free vibration of prestressed elliptic cylinders. 

Quite recently, Sewall et al. (refs. 6 and 7) applied the Rayleigh-Ritz technique to the 
vibration analysis of elliptic cylinders having simply supported, free, and clamped edges. 
Culberson and Boyd's (ref. 8) vibration analysis of simply supported oval cylinders 
involved a Fourier approach similar to that employed by Romano and Kempner (refs. 9 
and 10) in the static stress analysis problem. 

The present study was motivated by the approximate character of the solutions 
obtained, the lengthy algebraic manipulations involved, and convergence difficulties 
reported in the publications just cited. The purpose of this paper is to present a simple 
and accurate "mixed" finite -difference scheme for the determination of the natural fre- 
quencies and mode shapes of simply supported (longitudinal displacement is unrestrained), 
isotropic, closed arbitrary cylindrical shells with constant or variable wall thickness. 

The term "mixed" refers to the fact that both stress resultants and displacements are 
chosen as primary variables. The accuracy of the scheme is demonstrated by means of 
a number of numerical examples. 

The analytical formulation is based on the linear Sander s-Budiansky first- 
approximation shell theory. (See ref. 11.) A Fourier approach is used to separate the 
variables, and the governing equations are reduced to eight first-order ordinary differ- 
ential equations in the circumferential coordinate. The reduction to first-order differ- 
ential equations is similar to that used by Goldberg (ref. 12), Kalnins (ref. 13), and others 
in conjunction with numerical integration techniques. In contrast to the derivations in 
references 12 and 13, however, the governing equations in the present study are arranged 
to yield a symmetric coefficient matrix for the differential equations. For cylindrical 
shells with variable geometric or elastic characteristics, the symmetry of the governing 
equations is a unique feature of this formulation and is particularly convenient when a 
finite-difference method is used for the solution, since the resulting finite-difference 
equations are also symmetric. 

The finite -difference scheme presented herein is a modification of the method used 
in references 11 and 14 for the static analysis of the same class of shells. It is based 
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on the use of two interlacing grids for the different fundamental unknowns in contrast to 
other methods, including references 12 and 13, in which all unknowns are specified at the 
same set of points. References 12 and 13 use numerical integration (and multisegment 
method). The first-order equation formulation simplifies the finite -difference discreti- 
zation, and the use of interlacing grids results in reducing both the local discretization 
error and the bandwidth of the finite -difference field equations. In addition, the use of 
the interlacing grids results in a monotonic convergence for both the frequencies and the 
mode shapes with both reduction in local discretization error (through the use of higher 
order difference-quotient expressions) and mesh refinement. The cited advantages, 
coupled with the symmetric nature of the equations, lead to considerable improvement in 
the computational efficiency of the finite-difference scheme. 

The formulation presented herein is limited to cases where the fundamental 
unknowns are separable with respect to the independent variables. This is the case for 
noncircular cylindrical shells having simply supported curved edges and mechanical 
characteristics independent of the longitudinal coordinate. The formulation can be 
applied to other boundary conditions through the use of Kantorovich's method in conjunc- 
tion with the Hellinger-Reissner mixed variational principle (ref. 15). Although the pre- 
sentation is focused on noncircular cylinders, the finite-difference scheme can be extended 
to other classes of structures. 


SYMBOLS 


a,b 

Co 


c 


semimajor and semiminor axes of oval (or elliptic) profiles, respectively 

reference extensional rigidity of shell, chosen to be average extensional 

Eh 0 

rigidity over length Lg, 2 


length ratio, Lg/^i 



E Young’s modulus of elasticity 

f circular frequency of vibration of shell, hertz 

h local thickness of shell 
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reference thickness, chosen to be average thickness over length 


^o 



h v L 
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l 


M 1 > M 2’ M 12 


m 

N 


N lf N 2 ,N 


12 


N 


12 


n 


n 

Qj>Q2 

R 


t 

u,v,w 

V 2 

x l’ x 2’ x 3 

3=R 0 / L i 


axial and circumferential lengths of shell, respectively 
finite -difference interval in ^2 -direction 
moment stress resultants 
number of longitudinal half-waves 
exponent of order of local discretization error 
direct stress resultants 
modified (boundary) stress resultant 

number of finite -difference intervals in (^-direction in quarter of profile 

number of circumferential waves 

transverse shear stress resultants 

local radius of curvature of middle surface of shell 

reference radius of curvature, chosen to be radius of a circle with same 
circumference as that of noncircular profile 

time 

displacement components of middle surface in coordinate directions 
effective (boundary) transverse shear stress resultant (see eq. (3b)) 
curvilinear coordinates (see fig. 1) 
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0 - h o/ R o 

6 amplitude of thickness variation, measured from average thickness over 

length L 2 

y measure of oval eccentricity 

X = h/h 0 

v Poisson's ratio 

oj angular frequency, radians per second 

2 rotation components of middle surface 

density of shell material 
curvature parameter, R^R 

angle made by the normal at any point to the minor axis (see fig. 1) 

X 

dimensionless coordinates on shell middle surface, % a = for < 2 = 1,2 
row matrix 
column matrix 

rectangular or square matrix 
{^ml’I^m} vectors of fundamental unknowns 

Superscript T indicates transpose of a matrix. 

MATHEMATICAL FORMULATION 
Shell Geometry 

The geometric characteristics of the shell are shown in figure 1 and are defined as 
follows: L^ and Lg are axial and circumferential lengths of the middle surface of the 
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Figure 1.- Shell geometry. 


shell; h and h 0 are the local and reference shell thicknesses (h Q is chosen herein 
to be the average thickness of the shell over the length L 2 ); a and b are the semi- 
major and semiminor axes of the cross section. 

In the present study, both oval and elliptic profiles are considered. The radius of 
curvature of the oval, in the form found in reference 9, is 

R ( X2 ) = “ 71 m 

1 + y cos 4 it — 

L 2 

where R Q is the radius of a circle having the same circumference as the oval and y 
is an eccentricity parameter, which is a measure of the noncircularity of the cross 
section, 0 i jyj < 1. 

The radius of curvature of the elliptic profile is given by 



where 9 is the angle made by the normal at any point to the minor axis (fig. 1). 

To cast the problem in nondimensional form, two dimensionless coordinates ^ 
and ^2 are introduced, where 

(a =M) 
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Also, the following dimensionless shell parameters are used: 


Hh)-t 0 


Governing Differential Equations 

The fundamental unknowns of the mixed formulation used herein are chosen to be 
the quantities that would appear in the statement of the boundary conditions along an edge 
parallel to the x^-axis. For the Sanders -Budiansky first-approximation theory, these 
quantities are the four generalized displacements u, v, w, and 0g and the four stress 
resultants Ng, Njg, Mg, and Vg, where 


M 


N 12 " N 12 


12 


2R 


(3a) 


V 2 = «2 + < 3b > 

Here Njg and M^g are the modified (symmetric) stress resultants used in the 
Sanders -Budiansky theory. (See ref. 16.) 

If the different shell stress resultants and displacements are expanded in a Fourier 
series in the longitudinal coordinate so that the simple -support boundary conditions along 
the curved edges are satisfied, then 
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vw^ 2 N 2 M 2 V 2 ] = £ [v m w m 0 2m N 2m M 2m V 2m ] sin nur^ e lwt 

m=l 


( 4 ) 


and 


uN 12 Q 


il= ii> N 

m=l 


12m Q lm 


cos e 


iwt 


where o> is the angular frequency in radians/ second. 

The governing equations for arbitrary cylindrical shells then reduce to eight first- 
order ordinary differential equations for each Fourier harmonic. For the mth Fourier 
harmonic, these equations can be arranged to yield a symmetric coefficient matrix as 
follows (ref. 14): 


A ; 0 

f H 1 



0 ! B 

X^mJ 



( 5 ) 


where 



W T - l^mVC,, M 2m /C 0 u m W m 


M T * |> *2m L l N 12m L l/ C c V 2m L l/ C o 

The [A] and [B] are 4 by 4 symmetric submatrices, pjj is a diagonal submatrix, 
and [m^ and [m 2 ] are diagonal mass matrices. The coefficients of these matrices 

are given in appendix A. 

Equations (5) are free of the derivatives of elastic and geometric characteristics 
of the shell. Therefore, in cases where these quantities vary, no numerical approxima- 
tion is needed, and, consequently, the overall accuracy of the solution should increase. 
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METHOD OF SOLUTION 



Finite -Difference Discretization 


For the application of the finite -difference method, it is convenient to divide equa- 
tions (5) into two groups and rewrite them in the following form: 

+ [ : i]i 


(6a) 

-w{4 - [h 

{dH m } = ^ 2 |Tn 2 ]^ m } 

(6b) 


The basic idea of the modified finite -difference scheme used in the present study is 
to define the first derivatives of each of the fundamental unknowns at points lying midway 
between the points of definition of the same unknowns. This can be accomplished by using 
two sets of interlacing grids for the two groups of fundamental unknowns <H m > and 

KY 


The modified schemes to be used in conjunction with equations (6) are obtained by 
defining and |dH m j at points lying between the points of definition of |H m | 

{d^m} as shown in figure 2. Consequently, the two groups of equations (eqs. (6a) 


and 


and (6b)) are satisfied at different sets of points. The use of the above-mentioned inter- 
lacing grids makes it possible to avoid averaging in the difference -quotient expressions 
of the first derivative (ref. 17). Such averaging can lead to a significant reduction in 
accuracy as well as to nonmonotonic convergence. This will be discussed in the numer- 
ical examples. 





! 


° Control points for jH^J and {d^ m | 
x Control points for J^J and {dH m J 

Figure 2.- Control points for fundamental unknowns 
in modified scheme. 
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The finite -difference equations which simulate the governing equations (eqs. (6)) 
are obtained by replacing the first derivatives in these equations by their appropriate 
difference-quotient expressions. For evenly spaced intervals, the difference -quotient 
expressions with orders of local discretization error ranging between 0(1 2) and 
where l is the finite -difference interval, are given in appendix B for both the modified 
and the conventional schemes. Herein "by conventional schemes" means the schemes 
with all the fundamental unknowns and their derivatives defined at the same sets of points. 

As an illustration, the ^-modified finite-difference equations at generic interior 
points i and i + 1 are given by 





i+l 



(7a) 


and 


- IK.}, * !K} lt2 - sKL) - “ 2 hL KL p»> 


The resulting finite -difference equations can be represented in the following com- 
pact form: 


[K]{Zf = -w 2 [m]|Z} (8) 

where [K] and [in] contain the "generalized" stiffness and mass distribution and 
jzf is the vector of unknowns composed of the subvectors |H m | and (j = i 

and i ± 1 in the conventional and the modified schemes, respectively) at the various 
finite -difference stations. 

Equations (8) are banded and for this formulation will be symmetric in both the con- 
ventional and the modified schemes. The total bandwidth depends on the order of the 
approximation used and equals 8(N + 1) + 1 and 8(N - 1) + 1 in the conventional and 
the modified schemes (ref. 14), respectively, where N is the exponent of the local dis- 
cretization error (appendix B). 

Eigenvalue Extraction Technique 

A variant of the inverse -power method with shifts similar to that presented in ref- 
erence 18 has been used for the determination of the natural frequencies and mode shapes. 
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The method is described in detail in reference 18. In the present study, advantage was 
taken of the symmetric banded form of the matrix [K], and a direct Gaussian elimination 
procedure was used for each iteration to evaluate the new trial vector {Z}. 

NUMERICAL STUDIES 
Accuracy and Convergence Studies 

Since the efficiency, among other factors, of the modified finite -difference scheme 
presented in this paper can best be assessed by its accuracy in the most unfavorable cases 
within the class of applicable problems and by the simplicity of its application, the scheme 
was applied to a large number of oval and elliptic cylindrical shells with increasing degree 
of noncircularity, having both uniform and variable thicknesses. The effects of reduction 
in local discretization error and of mesh refinement on the accuracy and rate of conver- 
gence of the natural frequencies and mode shapes from the modified scheme were studied 
and compared with those from the conventional schemes. 

The class of shells considered in the present study included the simply supported 
elliptic shells analyzed by Sewall, Thompson, and Pusey (ref. 6) and the oval shells 
studied by Culberson and Boyd (ref. 8). 

In addition, shells with variable thickness were studied, in which case the thickness 
variation was assumed to be of the form 

h = h Q (l + 6 cos 47 r (9) 

where 6 represents the amplitude of thickness variation measured from the average 
thickness over the length L 2 - All the shells analyzed had doubly symmetric profiles, 
and consequently, only one -quarter of the shell circumference was considered with the 
boundary conditions at the ends taken to be the symmetric or skew -symmetric conditions. 
The four different combinations of symmetry and skew symmetry of the mode shapes, 
along the axes of symmetry of the profile, were studied. 

The local discretization error was reduced successively from O(fi) to O(f*0), 
and the number of finite -difference intervals in the shell quarter n was increased from 
10 to 40. The reduction of local discretization error up to 0(1 was feasible because 
of the extreme simplicity of the form of the governing differential equations (eqs. (6)). 

The numerical discretization was further simplified by the absence of boundary discreti- 
zation error in the closed shells considered in the present study. The numerical results 
presented herein pertain to the minimum frequencies and some of the associated mode 
shapes. 
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Before conducting the convergence studies, the minimum natural frequencies for 
one and two longitudinal half-waves were obtained by the finite-difference method and 
compared with the analytical results from references 6 and 8 for 6061 aluminum cylinders 
having the following characteristics: = 61.0 cm (24 in.), = 191.52 cm (75.4 in.), 

h = 0.0813 cm (0.032 in.), E = 68.95 GN/m2 (io 7 psi), v = 0.3, p s = 2768 kg/m 3 
(2.588 x 10"4 lb-secVin^). These characteristics were used throughout the numerical 
studies presented in this paper. Both the conventional and the modified solutions con- 
verged to the same frequency, and the converged results are summarized in table 1, 
where it is seen that the agreement between the finite-difference results and those of 
previous investigations is excellent. This agreement is particularly significant, since 
the finite -difference scheme presented herein is very simple to apply. This simplicity 
applies to the formulation of the problem, discretization, and computer implementation. 
Moreover, the difficulties with convergence as the degree of noncircularity increases 
(ref. 6, p. 14) or for higher Fourier harmonics (ref. 8, p. 1479) were, as is shown sub- 
sequently, totally nonexistent in the modified scheme presented herein. 

To elaborate on this last point, the convergence studies of the minimum natural 
frequencies obtained for two shells with high degree of noncircularity are summarized 
in table 2. In addition, figures 3 to 5 provide some indication of the accuracy and rate 
of convergence of both the modified and the conventional schemes as applied to a repre- 
sentative problem of an oval shell with y = 0.70 for three different modes. In fig- 
ures 3 and 4, the rate of convergence with both reduction in local discretization error 
and mesh refinement is shown for the three modes. Figure 5 shows the mode shapes 
associated with the minimum frequency in the three cases. 

The results of these numerical studies can be summarized as follows: 

The natural frequencies and mode shapes obtained by the modified scheme were, 
in general, more accurate than those obtained by the corresponding conventional schemes 
having the same order of local discretization error and the same number of finite- 
difference intervals. In addition, the bandwidth of the resulting finite -difference equa- 
tions is less in the modified scheme than in the corresponding conventional scheme. The 
increase in accuracy of the predictions of the modified scheme is mainly attributed to the 
use of "unaveraged" difference -quotient expressions. An indication of this accuracy is 
provided by the smaller error term in the first term of the remainder in the modified 
difference -quotient expressions (appendix B). 

In general, the difference between the predictions of the modified and the conven- 
tional schemes decreases with the reduction of the local discretization error. 

Although the convergence of the natural frequencies obtained by the modified scheme 
with both mesh refinement and reduction in local discretization error was monotonic, the 
convergence of the conventional schemes was, in some cases, nonmonotonic. (See figs. 3 
and 4.) 
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10 15 20 25 30 

n, Number of pivotal intervals 


(b) m = 3- 



N, Order of local discretization error 


(c) in = 6. 

Figure 3.- Convergence of minimum frequency 
with reduction in local discretization 
error. Oval shell with y = O.'JO, 



n. Number of pivotal intervals 

(c) m = 6. 

Figure 4.- Convergence of minimum fre 
quency with mesh refinement. Oval 
shell with y = 0.70. 







1.0 



Figure 5.- Mode shapes associated with minimum frequency for 
oval shells with 7 = O.7O. 


For all the homogeneous shells analyzed, the frequencies obtained by use of the 
Z^-modified scheme were identical with those obtained by use of the -conventional 
scheme with twice as many finite -difference intervals. This is not a coincidence, since 
the finite -difference equations for the Z^ -modified scheme are identical with those of the 
corresponding conventional scheme with twice as many finite -difference intervals. 

If the local discretization error is successively reduced from 0{l%) to 
the gain in accuracy becomes less appreciable for each order of reduction, the more so 
the finer the mesh used. This is especially true for the modified schemes, where the 
gain in accuracy obtained by reducing the local discretization error beyond 0(Z®) is 
almost insignificant. Since the penalty of increasing the bandwidth for each order of 
reduction in the discretization error remains constant, the use of modified schemes with 
discretization error of order higher than z® is not recommended unless the algorithm 
used for handling the algebraic equations does not depend on the bandwidth. 

Parametric Studies 

After the reliability of the results obtained by the modified schemes had been 
established, a limited number of parametric studies were made to provide some insight 
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into the effects of variations in the degree of noncircularity and local thickness on the 
minimum frequencies of the same group of oval and elliptic cylinders considered in the 
section entitled "Accuracy and Convergence Studies." 

The ratios a/b (major to minor axis) and h max /h m i n (maximum to minimum 
shell thickness) were taken as the measures of noncircularity and thickness variation, 
respectively. The thickness variation was assumed to be given by equation (9). Such a 
thickness variation for a positive 5 amounts to increasing the thickness in the shallow 
regions (ends of the minor axis) and decreasing it in the deep regions of the surface 

(ends of the major axis). The ratio h m ax/hmin is given by j-- ft . The ratios a/b 

and h maX| /h m i n were varied between 1.0 and 2.0. 

The effect of variation of a/b on the minimum frequencies of vibration is given 
in figure 6. It is shown that the minimum frequency of vibration decreases with the 
increase in a/b, that is, with the increase in the degree of noncircularity. This effect 
is more pronounced in oval shells than in elliptic shells having the same a/b ratio; 
this is especially true at higher values of a/b (a/b g 1.5). 



f = Minimum frequency of the noncircular cylinder 
f = Corresponding frequency of circular cylinder 

Figure 6.- Effect of degree of noncircularity on minimum frequencies of 
simply supported cylinders for one and two longitudinal half-waves. 
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As an indication of the effect of variation in thickness on the minimum frequencies 
of vibration, table 3 gives the frequencies for both oval and elliptic cylinders having 
a/b = 2 for different values of h m ax/ h min- The results presented in this table show 
that the chosen thickness variation resulted in an increase in the minimum frequency 
with increase of the ratio of h max /h m i n . The increase in frequency was somewhat 
larger for oval shells than for elliptical shells. Neither, however, changes the minimum 
frequency by more than about 12 percent. 

CONCLUDING REMARKS 

In this paper a modified mixed finite -difference scheme has been presented for the 
free -vibration analysis of simply supported closed arbitrary cylindrical shells with uni- 
form and variable wall thickness. The scheme is based on the linear Sanders-Budiansky 
theory, and the governing equations consist of eight first-order ordinary differential 
equations in the circumferential coordinate. In the finite -difference discretization, two 
interlacing grids are used for the different fundamental unknowns in such a way as to 
reduce both the local discretization error and the bandwidth of the resulting finite- 
difference field equations. The inverse-power method is used to determine the eigen- 
values and eigenvectors. 

Numerical studies were made of the effects of reducing the local discretization 
error (from order Oft 2 ) to order Oft 1 ®) where l is the finite -difference interval) 
and of mesh refinement on the accuracy and convergence of solutions obtained by the 
modified scheme presented herein, as well as by the conventional scheme (with all the 
fundamental unknowns defined at the same set of points). In addition, parametric studies 
were made of the effects of variation of the degree of noncircularity and local thickness 
on the minimum frequencies of oval and elliptic cylinders. These parametric studies 
show that the reduction in the minimum frequency of vibration with the increase in the 
degree of noncircularity is more pronounced in oval shells than in elliptic shells having 
the same major to minor axis ratio. Also, a slight increase in the minimum frequency 
can be obtained by means of a redistribution of the shell thickness. 

The finite -difference scheme presented herein is shown to combine a number of 
advantages over other finite -difference schemes previously reported in the literature. 
These advantages include the simplicity of the form of the governing differential equa- 
tions, the absence of the derivatives of elastic and geometric characteristics in these 
equations, and the symmetry of their coefficient matrix. Although the first two advan- 
tages result from the use of first-order equation formulation and are independent of the 
method of discretization, the symmetry of the coefficient matrix can only be efficiently 
utilized in a finite -difference discretization. The cited advantages result in simplifica- 
tions of the finite -difference discretization and provide flexibility for improving the 
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accuracy of the scheme without undue complications. Moreover, the effort required in 
the computer implementation of the scheme (coding, debugging, and verification) can be 
expected to be significantly less than that required for other finite -difference schemes. 

The numerical studies demonstrate the high accuracy of the predictions of the mod- 
ified scheme, even when a small number of finite -difference intervals are used, and the 
monotonic character of the convergence of the natural frequencies with both reduction in 
local discretization error and mesh refinement. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 14, 1972. 
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APPENDIX A 


FORMULAS FOR COEFFICIENTS IN THE GOVERNING EQUATIONS 
The independent nonzero terms of the submatrices [a] and [b] in equations (5) 


are given as 


A ll = "x 


Ajg = -muttc 


A = 
A 14 5 


* -12c 

22 ' X 3 ^ 2 

A 2 4 = wn 2 7r 2 c 


A 33 = m 2 7T 2 cX(l - iP) 


A 44 = 


i 4 jr^cX 3 a 2 i3 2 (l - iP) 


B 13 = -m7rc 




B 22 = ' 


m 2 77 2 a 2 /3 2 X 3 c(l - v) 


6 1 + 


177Q!pX 2 (P ( 


B 23 = / 922 


12 1 + 
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APPENDIX A - Concluded 


B 24 " ' c 


B 33 " 


2c 


/ ft 2 > 2 n 2 

x ( i-v)(i + 4^ 


is a diagonal submatrix given by 


M- 


-l 


-l 


[0] is a 4 by 4 null submatrix, and jmj and jmj are 4 by 4 diagonal mass sub - 
matrices given by 


ro 



0 


1 


1 



0 


where 


K = - 


Ps Xh o cL l 


and p c 


is the mass density. 
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APPENDIX B 


DIFFERENCE -QUOTIENT EXPRESSIONS USED IN PRESENT STUDY 

For evenly spaced intervals, the difference -quotient expressions at interior points 
for the conventional and the modified schemes can be written in the following compact 
form (fig. 7): 



a F- + 
±q i±q 


E 


R 


(Bl) 


where q = 0, 1,2,3,. . . for the conventional scheme, and 

dF li = f I a ±(l +4 )/2 V + E R <B2) 

q 

where q = 1,3, 5, 7,. . . for the modified scheme. In equations (Bl) and (B2), F denotes 
any of the fundamental unknowns |H m j> or j-; a r denotes weighting factors with 

a = -a r and aQ = 0; l is the finite -difference interval; Ej^ is the first term in the 
error series and is given by (ref. 19) 


E R =K R* 


N d N+1 F 


dx 


N+l 


,N ; 


where Z is the order of approximation and Kj^ is a multiplier. The ± sign before 



° Control points for F and dF 

(a) Conventional schemes. 



o Control points for F 
x Control points for dF 

(b) Modified schemes. 


Figure 7-- Control points for F and dF in conventional and modified schemes. 
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APPENDIX B - Concluded 

a term denotes that the summation extends over the positive as well as the negative 
values of the index. In each term the sign of the two indices of a and F should be 
identical. 

The different values of a r and K, ;> for both the conventional and the modified 
schemes with local discretization errors up to 0(^10) are given in table 4. 
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TABLE 1. - COMPARISON OF MINIMUM FREQUENCIES OBTAINED 
BY FINITE -DIFFERENCE METHOD WITH THOSE 
OBTAINED BY PREVIOUS INVESTIGATORS 


m 

h 

a/b 

f, Hz, for elliptic profile 

f, Hz, for oval profile 

Reference 6 

Finite 

differences 

Reference 8 

Finite 

differences 

1 

7 

1.0 






7 

1.176 






7 

1.538 

138.5 

138.6 

129.6 

129.9 


6 

2.5 

106.8 

106.8 

— 

— 

2 

10 

1.0 

325.1 

325.7 

327.1 

325.7 


10 

1.176 

310.6 

310.7 

310.4 

309.3 


10 

1.538 

268.1 

268.2 

245.7 

245.4 


9 

2.5 

202.4 

202.3 


— 


a Skew-symmetric mode, all others are doubly symmetric mode. 
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TABLE 2.- CONVERGENCE OF MINIMUM FREQUENCIES OBTAINED 
BY DIFFERENT FINITE -DIFFERENCE SCHEMES - Concluded 


(b) Elliptic cylinder with 1.538, m = 3 



Frequency, Hz, for - 

n 

I 2 

I 4 

Z 6 

z 8 

z 10 

Modified schemes 

10 

380.9 

389.5 

392.2 

393.4 

393.9 

15 

389.4 

393.5 

394.3 

394.5 

394.5 

20 

391.8 

394.2 

394.5 

a 394.6 

a 394.6 

25 

392.8 

394.4 

394.5 



30 

393.3 

394.5 

394.5 

* 



Conventional schemes 

10 

563.8 

406.5 

375.9 

374.5 

380.6 

15 

374.9 

386.6 

392.1 

393.7 

394.3 

20 

380.9 

392.2 

394.1 

394.5 

394.5 

25 

386.8 

393.6 

394.4 

394.5 

a 394.6 

30 

389.4 

394.1 

394.5 

394.5 



a Converged solution. 




TABLE 3.- EFFECT OF THICKNESS VARIATION ON MINIMUM 
FREQUENCIES IN OVAL AND ELLIPTIC CYLINDERS 
WITH ^ = 2, m = 2 


^max/^min 

Minimum frequency, Hz 

Elliptic profile 

Oval profile 

1.0 

230.1 

107.3 

1.25 

237.6 

112.0 

1.50 

242.9 

115.5 

1.75 

246.9 

118.4 

2.0 

250.0 

120.7 
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TABLE 4.- COEFFICIENTS IN THE FINITE -DIFFERENCE -QUOTIENT EXPRESSIONS 
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